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We compute the angular power spectrum Ci of the BATSE 3B cata- 
log, and find no evidence for clustering on any scale. These constraints 
bridge the entire range from small scales, probing source clustering and 
repetition, to large scales constraining possible Galactic anisotropies, 
or those from nearby cosmological large scale structures. 



Q\ I. INTRODUCTION 

ON 

X) 

The observed angular distribution of 7-ray bursts (GRBs) is isotropic, while 
["T 1 their brightness distribution shows a reduced number of faint events. These 
observations favor a cosmological burst origin. Clustering of bursts could 
be evidence of actual clustering of sources or of repeated emission. Repeti- 
tion would call into question the viability of many cosmological burst models. 
^} Anisotropies manifest themselves on different angular scales and with differ- 
ent magnitudes. Galactic features cause large-scale distortions, while true 
repetition would affect small scales. For large-scale signatures, we search for 
excesses of sources towards some direction or a concentration towards some 
^sO plane in the sky, i.e., we seek a dipole- or quadrupole moment. It is now 
ON common practice to apply both coordinate-free and galactic tests (1). Dipole- 
and quadrupole measures were sufficient when sample sizes were small. Now 
O an extension of moment methods to higher orders is needed. Low order multi- 
■ poles are not sensitive to instrumental smearing, but higher harmonics are. If 
^~ associated with galaxies, we expect clustering on very small scales. If bursts 
repeat, we expect clustering at 6=0. Both effects are diluted by localization 
uncertainties, and angular power is transferred from small (or zero) angular 
scales to a scale given by the detector response. One tool for the analysis of 
source clustering is the two-point correlation function (2), which is related to 
the power spectrum through a Fourier transform (4). 
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II. METHOD 



We model the GRB distribution as a 2D stochastic point process n(r) = 
S(r, r 8 -) with intensity (average point density per steradian) A(r). Here S 
denotes the 2D Dirac delta function, and the unit vectors r 8 - correspond to 
the various GRB positions. If we had detected a nearly infinite number of 
bursts, then the function A(r) would be known with great accuracy, and the 
only source of errors when computing its power spectrum would be cosmic 
variance. Since in practice we have only a finite number of bursts (1122 for 
3B), our estimates of A include shot noise. A Poisson process satisfies the 
expectation value equations 

(n(r)) = A(r), (1) 

and 

(n(r)n(r')) = A(r)A(r') + S(r, r')A(r). (2) 

Here A is itself a random field, A(r) = n(r)[l + A(r)], where the underlying 
density fluctuations A are modeled as a Gaussian random field. The function 
n, which we will refer to as the exposure function, is the number of bursts 
per steradian expected a priori, not the number density actually observed. In 
other words, n(r) is proportional to the exposure time in the sky direction r. 
We assume that (A(r)) = and that the statistical properties of the field A 
are isotropic, which means that if we expand it in spherical harmonics as 

oo I 

A ( r ) =J2 J2 a ^Yi m {r), (3) 

1=0 m = -l 

then 

{aim^l'm') = dll'dmm'Cl, (4) 

where the coefficients C'i are known as the angular power spectrum. There are 
thus two separate random steps involved in generating n: first the generation 
of the smooth field A, then the Poissonian distribution of points. 

Given the field n(r), we wish to estimate the coefficients ai m . We define 
them as 

f n ( r ) 

atm = Y lm (r)——dCl-6 l0 6 m0 \/4:ir. (5) 



We now compute the statistical properties of these estimates. By substitution 
we obtain 

{aim) = J Yim{r)di} - SioSmoV^TT = 0, (6) 

i.e., the expectation values vanish. Since the expectation values of the true 
coefficients ai m vanish as well, this means that our estimates are unbiased. 



Using the expressions above, we find that the correlation between two multi- 
pole estimates is 



{ai m aii m i)= I I Ytm(r)Yfim>(r') 
which reduces to 



(A(r)A(f')> + -U(r,f') 



dndn', (7) 



I- - \ x x r< , / Y z m (r)Yz, m ,(r) 
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Defining the quantities 

Clm = aim ~ blm, (9) 

we find that they are unbiased estimates if we choose the bias correction as 

hm = [ (10) 

If n is constant, then the bias correction becomes simply bi m = l/«, indepen- 
dent of I and m. The C'i m are thus good estimates of C'i for each m- value 
separately. To reduce error bars, we estimate power by averaging the Ci m - 
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Ct = 21+1 ^ ^ tm ' ^ 

m — — i 

Defining b to be the average of the bias corrections bi m , we find that b is in 
fact independent of i, and obtain 

i A , i r drt 

21+1 4tt J n(r) 

m — — t 

i.e., b is just the spherical average of 1/n. 

It is straightforward to include the effects of position errors in the formalism, 
which is described in a more detailed ApJ version of this paper (5). We model 
the BATSE beam function as a Fisher function 

exp [<7 _2 r • r'l , % 

B(r-r') = A \ . ur ^ 13 

v ; 4ira 2 smhfo- 2 ] v ; 

characterized by a location error a. This is a spherical version of the Gaussian 
distribution, and reduces to 

exp f-i£ 

B(co S e)^ 1 2 J (14) 

when cr < 1 radian « 60° . The Fisher function has the advantage that it is 
correctly normalized (its integral over the sphere is unity) for arbitrarily large 



angles a, which is not the case for the plane Gaussian. In addition to statistical 
position errors we include (in quadrature) a 1.6° systematic uncertainty. This 
value is significantly lower than the 4° of earlier catalogs, allowing us to extend 
spherical harmonic analysis to I ~ 60 before localization uncertainties wash 
out intrinsic angular power. 



The power spectrum C'i of the the 3B data (3) is shown in Figure 1. There is 
no evidence of deviations from isotropy on any angular scale. If the gamma-ray 
bursts are completely uncorrelated, the points should scatter symmetrically 
around zero, with about 68% in the shaded region. Since all power is by 
definition positive, the presence of any type of clustering would shift the dis- 
tribution upwards, leading to a positive excess. A monopole Co /Ait = 0.0001 
corresponds to a fluctuation of VO.0001 = 1% in the average burst density. 
Likewise, [dj Att] 1 / 2 can be interpreted as the density fluctuation on the an- 
gular scale 9 « 60°/£. The size of the error bars (the height of the shaded 
region) is readily understood. For I = 0, all N = 1122 bursts carry equal 
weight, so apart from a factor of y/2, the shot noise gives just the familiar 
Poisson variance l/N. As I increases, the error bars become smaller since 
(2£ + 1) independent modes are being averaged. Since the weighting scheme 
loosely speaking only obtains information on C'i from bursts better localized 
than 60°/£, the l/N shot noise finally causes the error bars to grow with I 
again, since the effective number of bursts N decreases. 

Although the angular power spectrum C'i provides a useful measure of the 
amount of clustering on different angular scales, it does not contain any in- 
formation about the relative phases of the different multipoles at m . The loss 
of phase-information means that although the power spectrum may tell us 
that there is extra power on some scale, it does not tell us anything about 
where in the sky this power is coming from. Fortunately, this type of infor- 
mation is easy to extract through definition of a map xt(i), the multipole map 
corresponding to multipole i, as the sky map 



The multipole information that our SHA extracts from the data, as plotted 
in Figure 1, places sharp quantitative limits on repetition. Suppose that a 
fraction / of all observed bursts can be labeled as repeaters that are observed 
to burst v times each. Application of an SHA-based technique to test this 
two-parameter family of models against the BATSE 3B data shows that all 
models with (y — 1)/ > 0.05 are ruled out at 99% confidence (6), as compared 
to the best previous 99% limit (y— 1)/ > 0.27. Thus even a cluster of 6 events 
from a single source would have caused excess power above that present in 
the 3B catalog. 



III. RESULTS AND DISCUSSION 
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FIG. 1. The shot-noise corrected angular power spectrum of 3B (solid squares). 
The shaded region shows the la error bars. Any type of clustering would drive the 
measured points upward. The double-shaded region shows what the errors would 
be without localization uncertainties. 

In summary, multipole expansion of the projected distribution of GRBs does 
not show evidence for clustering on any angular scale. This argues against 
the recurrence of a substantial fraction of burst sources (6) and against any 
source population with strong intrinsic anisotropies. The remarkable degree of 
isotropy of GRBs severely constrains any burst model that invokes traditional 
geometric features of the Milky Way (disk, bulge, or halo). 
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